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PREFACE 



This report constitutes a pre-presentation publication 

# * V* 

of a paper and results to be presented at the GAMM Workshop 
on ’'Numerical Methods for the Computation of Inviscid Trans- 
onic Flow with Shock Waves". This conference will be held 
at the Aeronautical Research Institute (FFA), Stockholm, 
Sweden on September 18-19, 1979. 

This report includes the conference paper, all results 
obtained for the various workshop test cases, and approximate 
iso-mach line plots for the supercritical cases. 

This work was primarily supported by NASA Grant NSG1174 
and partially supported by the Texas Engineering Experiment 
Station. 
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I . Introduction 


This paper will briefly discuss some of the results ob- 1 
tained in the process of solving the test problems for the 
GAMM Workshop on "Numerical Methods for the Computation of 
Inviscid Transonic Flow with Shock Waves" with the TRANDES 
program. Briefly, this method, 1 " 5 utilizes the full, invis- 
cid, perturbation-potential flow equation in a Cartesian grid 
system that is stretched to infinity. This equation is rep- 
resented by a non-conservative system of finite difference 
equations that includes at supersonic points a rotated differ- 
ence scheme and is solved by column relaxation. The solution 
usuaJly starts from a zero perturbation potential on a very 
coarse grid (typically 13 x 7) followed by several grid halv- 
ings until a final solution is obtained on a fine grid (97x49). 
Occasionally, for cases having high local Mach numbers, the 
solution must be started on the coarse grid (25 x 13). Since 
the airfoil does not coincide with the grid points, the sur- 
face boundary conditions are represented as two-term Taylor 
series about dummy points inside the airfoil. On the outer 
boundaries, the exact infinity conditions are used. This 
method can, if desired, include the effects of weak viscous 
interaction or be used in the design mode. 2-6 

All of the results presented at the workshop and in this 
paper were obtained at the rate of 10 4 pts/sec on an Amdahl 
470/V6 using a FORTG compiler and single precision arithmetic 
(less than 7 significant digits). A typical run took 4-8 min- 
utes, although good engineering results were obtained in one 
minute on the medium grid; and convergence was obtained on the 
meldium grid to at least a maximum cyclic perturbation change 
of less than 5E-5. On the fine grid, the A^max was usually 
larger due to significant digit error in the far-field. It 
should be noted that the workshop cases were run 1. 5-5.0 times 
longer than usual due to the availability of computer time 
from TAMU, and that the TRANDES program has never been opti- 
mized for time. 

In order to maintain the goal of common discretization, 
the grid stretchings were setup so that the number of airfoil 
points, number of wake points, Ax and Ay at the trailing edge, 
and the location of the last finite vertical grid column 
matched the suggested grids as closely as possible. For cases 
involving large supersonic zones, the y-grid was extended so 
that the last finite horizontal grid line was subsonic. Other- 
wise, the rotated difference scheme might have us^d undefined 
values . 
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for the test cases appear to show correct trends. (I.E. For 
C. = 0.0 — M = 0.72, CD = .0004; M = 0.8, CD = 0.0100; M = 

0.85, CD = 0.0381; M a = 0.95, CD = 0.0989) 

b. Bump in Channel 

Thi% problem was solved by treating the bump as a sym- 
metrical airfoil in a solid wall wind-tunnel. The solution 
used 41 points on the upper surface, 16 in the wake, 21 pts 
vertically from the centerline to the wall, AX te = 0.0251, 

AY. = 0.03, and <f> = 0 on the channel walls. For the up- 

y 

stream/downstream infinity conditions two approaches were 
tried. The first used the asymptotic form derived by Murman 8 
while the second had <j> =0 imposed and allowed the solution 

to float. The results x on the bump and between it and the 
wall were identical. However, some slight differences were 
observed upstream and downstream, with the floating solution 
showing less blockage type influence. While the floating 
approach is easier to implement, it is more prone to signif- 
icant digit errors due to the magnitude of the floating 
potential (0(1) instead of 0(0.1)). 

Finally, several other upstream/downstream boundary con- 
ditions were imposed, such as <<>=0; and all yielded essen- 
tially the same results. Since the channel case was coded 
rapidly for the workshop, this lack of sensitivity may be 
due to either coding errors or a special feature of the test 
case. 

c. RAE 2822 

This problem was solved using the same grid as for the 
NACA 0012 cases. During the solutions, it was discovered 
that the subcritical results were sensitive to the location 
of the first and last points on the airfoil, which are 
normally at 0.01c and 0.99c (X4 = 0.49). The initial results 
(solid line. Fig. 2) indicated a possible peak in C near 

the leading edge. To resolve this peak, the first points 
were moved to 0.005c and 0.995c (X4 = 0.495), and slightly 
different results were obtained as shown on Fig. 2. Inter- 
estingly, EPSS had to be increased to 1.0 for this second 
case due to the appearance of supersonic points on the 
coarser grinds. Obviously, the Cartesian grid placement 
sometimes affects the airfoil effective a; and it is probably 
best to Correlate results versus C T rather than a . 


d. CAST 7 

Again the basic NACA 0012 grid was used. Here, in both 
the subcritical and supercritical cases, oscillations were 
observed in C , at 0.03c and in C at 0.92c. Since the sur- 
face slopes w^re also oscillatory* 3 !^ these areas, it is be- 
lieved this behavior was due to the sparsity of the given 
coordinates and the use of spline fits to determine the comp- 
utational ordinates and slopes. Also, in both cases, the 


aft C i. 
Pi 

nation . 
(12.5°) 
trend. 


bucket turned around at the TE and approached stag- 

This behavior is reasonable due to the large TE angle 
In actuality, viscous effects would mitigate this 


e . KORN 1 


Th% basic NACA 0012 grid was also used for this case. 

For the design point, the resultant C distribution was close 
to the theoretical hodograph values aRd was almost shockless. 
Again some sensitivity to the grid size and X4 value was 
observed,,- as can be seen in Fig. 3. Also, a double precision 
calculation (16 digits) was made and essentially yielded the 
same results as with single precision (7 digits). However, 
for the 16 digit case the A^max on the fine grid was steadily 
decreasing instead of oscillatory and smaller (6E-5 vs. 2E-4). 

III. Conclusions 

Except for the NACA 0012, M m = 0.95 case, all test 
problems were solved straight-forwardly and appeared to be 
converged or close to convergence. The only difficulty was 
some sensitivity to grid placement, which is typical of the 
Cartesian formulation. Again, note that all results were 
obtained in single precision (less than 7 significant digits) 
on an Amdahl 470/V6. 
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Summary of Test Problem Results 
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Figure 1: Case A-I( i) , Airfoil Cp Distribution versus X/C 
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Figure 2: Case A-I(i), Airfoil C„ Distribution versus Y/C 










Figure 4: Case A^-I( ii) , Airfoil Cp Distribution versus Y/C 
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Figure 8 : Case A-II(ii) ,Airfiol C p Distribution versus Y/C 
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Figure 10: Case A-II(iii), Airfoil C p Distribution versus X/C, into wake region 






Figure 11: Case A-II(iii) ,Wake C Distribution versus X/C 
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Figure 12 : Case A-I I( iii) , Airfoil Cp Distribution versus Y/C 
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Figure 13: Case' A-II( iv) .Airfoil C n Distribution versus X/C 
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Figure 17: Case B, Upstream C n Distribution 
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Figure 18: Case B/'Bump" Distribution versus X/C 
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Figure 20: Case B, "Bump" C. Distribution versus Y/C 
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Figure 22: Case C-I , Airfoil C p Distribution versus Y/C 




Figure 23: Case C-II, Airfoil C Distribution versus X/C 





Figure 24: Case C- I I, Air foil C Distribution versus Y/C 
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Figure 25: Case D-I, Airfoil C Distribution versus $/C 
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Figure 27: Case D-II , Airfoil C„ Distribution versus X/C 
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Figure 29: Case G-I, Airfoil C Distribution versus X/C 
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